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Abstract 

o . 

fN| , Integrable discretisations for a class of coupled nonlinear Schrodinger (NLS) type of equa- 

tions are presented. The class corresponds to a Lax operator with entries in a Grassmann 
algebra. Elementary Darboux transformations are constructed. As a result, Grassmann 
generalisations of the Toda lattice and the NLS dressing chain are obtained. The cani- 
ty- | patibility (Bianchi commutativity) of these Darboux transformations leads to integrable 
Grassmamm generalisations of the difference Toda and NLS equations. The resulting 
discrete systems will have Lax pairs provided by the set of two consistent Darboux trans- 
formations. 
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1 Introduction 



There has been a constant interest to noncommutative extensions of integrable equations over 
the last few decades [6, 41, 9, 11, 15, 17, 24, 36, 37, 40]. Among the well-known examples are 
the noncommutative analogues of KdV, NLS, sin-Gordon, the KP equation, the Hirota-Miwa 
equation, two-dimensional Toda lattice equation and AKNS hierarchy. 

Supersymmetric systems are particular examples of noncommutative extensions of inte- 
grable systems. They have attracted much attention because of their applications in physics. 
Perhaps the best known example of such equation is the Manin-Radul super-KdV equation 
[30]. It has been extensively studied by various authors [35, 28, 29]. In [28] soliton solutions 
for the Manin-Radul super-KdV equation were constructed by iterating Darboux transforma- 
tions, considering the cases of an even number and an odd number of iterations separately. 
In [ "] unified formulae for the solutions not depending on the parity of the number of itera- 
tions are obtained. This is based on an alternative approach to Darboux transformations using 
quasi-determinants [14]. 

Moreover, the connections between the supersymmetric soliton equations and the theory 
of Lie superalgebras were analyzed in [16, 34, 19], in order to obtain a super-analogue of the 
classical Drinfeld-Sokolov theory [12]. 

Darboux transformations play a crucial role in the theory of integrable systems [31]. They 
are gauge-like transformations which preserve the form of the associated linear problem (Lax 
representation) and are basic tools for obtaining the so-called soliton solutions [43, 44, 9, 8]. 
The theory of Darboux transformations was boosted by the dressing method [42]. 

Darboux transformations also play a role in constructing integrable discretizations of in- 
tegrable equations [3]: to a very wide variety of (continuous) integrable nonlinear PDE's the 
associated dressing chains can be interpreted as Backlund-Darboux transformations [26]. Fur- 
thermore, the compatibility (Bianchi commutativity) of these Backlund-Darboux transforma- 
tions leads to equations with two discrete independent variables (lattice equations). The Bianchi 
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commutativity for Backlund-Darboux transformations is also known as a principle for nonlinear 
superposition [ ]. The classifications of elementary Darboux transforms can be used as a tool 
to classify discrete systems related to a given Lax operator. These discrete systems will have 
Lax pairs provided by the set of two consistent Darboux transformations. The corresponding 
Backhand transformations will represent symmetries of the discrete (difference systems). 

Here we will study Darboux transformation for a class of the nonlinear Schrodinger (NLS) 
type of systems with Lax pairs belonging to Grassmann algebras. Such systems are closely- 
related to the supersymmetric version of the NLS equation with an osp(l|2)-invariant Lax pair 
studied by P. P. Kulish [22]: 

iut + Uxx ~ 2u ] uu - + iW x = 

i^ 4 + V xx - v)v$> + \{2v$! x + &u x ) = 0' 

Here u and \I/ are smooth functions taking values in an infinite-dimensional Grassmann algebra 
S = So©Si- The variables u are called commuting (Bosonic) variables: u\u 2 = U2U1, ui,u 2 G So 
while the variables \I/ are called anti-commuting (Fermionic) ones: ^1^2 = — ^2^1, ^1, ^2 G Si- 
Here also f stays for a Hermitian conjugation in S- 

The paper is organised as follows: in Section 2 we present some basic facts of Grassmann 
algebras, the Lax representation and the general form of the partial differential equations. In 
Section 3 we describe the elementary Darboux transformations and the corresponding dressing 
chains. In Section 4 we present the discrete integrable systems obtained as a compatibility 
condition of elementary Darboux transformations from Section 3. 



2 Preliminaries: Grassmann algebras and Lax represen- 
tation 

Firstly, we present basic facts and properties of Grassmann algebras. Further details can be 
found e.g. in [4, 25]. 

Let S be a Z 2 -graded algebra over a field K of characteristics zero (such as C or Q). Thus, S 
as a linear space is a direct sum S = So © Si, such that SiSj Q Qi+j (mod 2). Those elements 
of S that belong either to So or to Si are called homogeneous, the ones from So are called even, 
while those in Si are called odd. 

By definition, the parity \a\ of an even homogeneous element a is and it is 1 for odd 
homogeneous elements. The parity of the product \ab\ of two homogeneous elements is a sum 
of their parities: \ab\ = \a\ + \b\. Grassmann commutativity means that ba = (-l)l a H fc la6 for 
any homogenous elements a and b. In particular, a\ = 0, for all a\ G Si and even elements 
commute with all elements of S- 

Consider a Lax operator of the form 

L = d x + U-Xh, (2.1) 
where the matrix U has entries in a Grassmann algebra: 

/0^2^\ /lfl 

U= -x C , h= - 
\2p <p J 2 \0 0-l 

Here p and q are even elements of the Grassmann algebra S, while £, x, <ft and ip are odd 
homogeneous ones and A G C is a spectral parameter (even variable). We will be using the 
natural grading Uij G Si+j (mod 2). 
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The zero curvature condition [L, A] = leads to the following system of equations: 

Qt = -Qxx + *Px( ~ i^Cx ~ 2q(ipx ~ <P() + 8q 2 p 
Pt = Pxx + o,-x - (pK x + 2p(ipx - 4>Q - 8p 2 q 

ipt = ipxx - Qx<P ~ 2q4>x + ^pqifj + tp(p( (2-3) 

Ct = Cxx - q x x - 2gx x . + 2pqC - 4>x( 

K t = -x xx + p x ( + 2p( x + 2pqn - (px( 
4>t = -<Pxx + Px^P + Zpipx + Zpq(j) + 1p(j)K 

The second Lax operator A is of the form: 

A — dt + Vo + XU - X 2 h, (2.4) 

where 

Upq-2ipx \ 

V = ad^ 1 U x + -Apq -2<f>( . (2.5) 

V 2(0C + ^x)/ 

The reduction ip = £ = (jy = w and p = q^ leads to a system which after a re-scaling and a 
point transformation t — > it, x — > ix leads to (1.1). It can be shown that the system (2.3) is a 
completely integrable Hamiltonian system. It has an infinite number of conservation laws. The 
first three constants of motion are of the form: 

N = JZc dx {Apq + <W> + xC}; (2.6) 

7 = I-oo dx {~ 2 P<? _ <M>* - H Cx - <Px1p - Xx(}] (2.7) 

^ = JZo dx { 2 Px<lx + 4>x^x + >*x(x + 2p 2 q 2 + pq{(pip + x£) 

-q(<p<p x + xx x ) - p(^ x + CCb)}- (2-8) 

In Physical applications N is known as a "total number of particles", 7 as a "total momentum" 
while IK is the Hamiltoniam of the system. 

3 Darboux transforms for the Lax operator (2.1) 

By a Darboux transformation we understand a map 

L^L X = MLM' 1 (3.1) 
where the Lax operator Li has the same form (2.1) but with an updated potential U\. 

( iP 1 2q 1 \ 

Lx = d x + U l - Xh, Ux = \-x x Ci • (3.2) 

V 2jh 0i / 

In (3.1) M is a matrix whose entries Mij G Si+j (mod 2) are rational functions of A and 
differentiable functions of x. It follows from (3.1) that 

M x + XJ X M -MU = 0. (3.3) 
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Obviously, a composition of Darboux transformations is again a Darboux transformation with 
more complicated rational dependence in A. 

We are interested in elementary Darboux transformations which cannot be decomposed 
further. Thus, we restrict ourselves by linear in A Darboux matrices: 

M = M + \M X . (3.4) 

The substitution of M in (3.3) results in: 

[h,Mx] = 0, (3.5) 

M M + ?7iM a - M X U + [h, M ] = 0, (3.6) 

M 0<x + UxMq - M U = 0. (3.7) 

Let us consider the simplest case of A- independent Darboux transformations (Mi = 0). From 
(3.6) it follows that Mq is a diagonal matrix. Then, equation (3.7) implies that Mq is a constant 
diagonal matrix, and that U\Mq = M Q U. The later is nothing but a Lie point symmetry 
transformation, which does not lead to non-trivial results. 

If Mi ^ 0, then it follows from (3.5) that Mi is a diagonal matrix: Mj = diag (a, /3, 7). 
Furthermore, the equation (3.6) implies that a, (3 and 7 are constants. We will describe the 
elementary Darboux transformations for the special case when the matrix Mi in (3.4) has rank 
one and a = 1, (3 = 0, 7 = 0. In this case it follows from (3.6) and (3.7) that M 0j22 and 
Mo j33 are constants. Further analysis shows that there are two essentially different cases: (1) 
M , 22 = and Mo, 33 = 1; (2) M , 22 = M , 33 = 1- 

For a sake of convenience, from now on, we will denote the matrix element Mq by F. 

Case 1: M 22 = 1 and M 33 = 0. The A-term of the compatibility condition (3.6) in this case 
gives the conditions: M i2 = ip, M 21 = — x 1; M 13 = p, M 31 = pi and M 23 = M 32 = 0. Therefore, 
the Darboux matrix M takes the form 



M— -xi 1 . (3. 




Then, the A-independent term in the compatibility condition (3.7) leads to the set of algebraic 
constrains: 

0i = -Piip, C = -q*u PiQ = !■ (3-9) 

and to the following system of equations: 

q x = -(-0Xi + 2F)q 

F x = 2 ( — — ) - *0ixi 

\Q Q-iJ 

^ = ^1-^ + — ^-1 (3.10) 

q-i 

>Ci x = —x + X!.r H x 2 

Now, if we introduce new variables v, <p and -0 together with the forward v± and backward u_i 
shifts of f 



? = e u , p = e"- 1 ! 



^ = ^Z 2 , X! = v^e"^ 2 , 
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using also the algebraic constrains (3.9) one can eliminate the function F and cast (3.10) into 
a system of 3 equations: 

v xx = 4 (e^ - e^" 1 ) + (p77_ a + tp^e^-^ 2 - (^i + ^)e^/ 2 
9x = ^e^-")/ 2 + v9_ ie (,; - v - l)/2 (3.11) 

The above system is an integrable noncommutative extension of the Toda chain: the reduction 
£ = r/ = leads to the standard Toda chain: 

Vxx = 4e vi ~ v -4e v - v -\ 

The system (3.11) also has a Lagrangian formulation with a Lagrangian: 

£, r/) = y dx ^ - 4e^- + 2(^_ 1 + (p^e^-^ 2 + cp Vx - <p x r^ . (3.12) 
Case 2: M 22 = 1, M 33 = 1. 

It follows from(3.6) that M 12 = ip, M 21 = -x i; M 13 = q, M 31 = p x and M 23 = M 32 = 0. 
Due to Abel's theorem, the Wronskian does not depend on x (since the potential U is a traceless 
matrix) and thus (F — p±q + ip*Ci) x — 0. Therefore, the function F can be determined up to a 
constant \i G So : 

F — piq — tpx l + pi. (3.13) 
As a result, the Darboux matrix M takes the form 

(^ + Piq- ip>c\ + X if) q\ 
-xi 10. (3.14) 

Pi l) 

Then, the condition (3.7) leads to the following algebraic equations: 

(l-T)C = (g*i), (l-T)<f>=( Pl i/j). (3.15) 

and to a set of 4 dressing chain equations: 

q x = 2q(tpxi - p x q - ji) + 2q 1 - (1 - Ty l (qxi)if) 
p x = -2p(if)-ix - pq-i - //) - 2p_i - (1 - T) _1 (pii/>)x 
ip x = i>i- 9(1 - T)~ l [piip) - {^ + P\<l)^ (3-16) 
= -x_i - p(l - T) _1 (gxi) + (/i + pg_i)x 

Here T is the shift operator induced by the Darboux transform (3.14). It acts on the potential 
as U\ = TU = AdMi^: Qi = Tq, p\ = Tp, ip\ = Tip <pi = Tcf), = T( and xi = Tx. The 
presence of the operator (1 — T) _1 in (3.16) leads to a non-local dressing chain. We should note 
that (3.16) can be rewritten into a local form in terms of the variables p, q, (f> and £, although 
it will lead to non-evolutionary dressing chain equations for the odd variables (f> and (. 

In the Bosonic limit (when all non-commuting variables vanish), (3.16) reduces to the stan- 
dard NLS dressing chain [3]: 

q x = -2q{p x q + jj) + 2q x p hx = 2p 1 (p 1 q + ji) - 2p. (3.17) 
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TS[V] = ST[V] 




Figure 1: Bianchi commuting diagram 

4 Darboux transforms and discretisation 

A standard approach for obtaining integrable difference equations is to construct a pair of 
Darboux transformations: discrete systems appear as consistency conditions of two Darboux 
matrices M and iV around a square [26, 3] (the Bianchi commutativity). 

Let us introduce lattice variables (k,m) (k,m e Z): generic even u&, m and odd variables 
Tk, m are defined on an integer lattice Z x Z (through the whole text we adopt the convention 
to denote by Greek letters the odd variables and with Latin letters the even ones). It is also 
useful to introduce the shift operators S and T as follows: the operator S shifts the first lattice 
variable while the operator T shifts the second one. For example 

Sqk,m = Qk+l,m, T(k,m = Tk,m+1, TS = ST. (4.1) 

Now, consider two Darboux transformations M(A) and N(X). On the space of fundamental 
solutions of L(X) (2.1) they act as follows: 

5[*(A)] = M(A)tt(A), T[tf (A)] = iV(A)*(A). (4.2) 

The compatibility of the transformations [S, T] = (Bianchi commutativity) implies 

S[N(X)]M(X) = T[M{\)]N{\), (4.3) 

and leads to a set of algebraic relations between U, S[U], T[U] and TS[U]. 

In this setting, Darboux transformations (4.2) can be considered as a discrete Lax pair 
associated with L(A), while the system of algebraic equations (4.3) can be considered as a 
system of difference equations for variables labeled by the point on the lattice. Hence, the 
system (4.3) is an integrable discretisation of the hierarchy of L(X) having a discrete Lax pair 
given by (4.2). Moreover, the differential equations for the Biicklund transformation, coming 
from the derivation of the Darboux matrices can be considered as symmetries of the difference 
system. 

Here we will describe the set of integrable discretisations obtained from imposing a consis- 
tency of two Darboux transformations, obtained in the previous Section. 

Case A: In this case, consider two Darboux matrices of the type described in case 2 in 
section 3: 

(A* + PioQ - V^io + A ip q\ ( n + poiq - i^oi + A ip q \ 

-xw 10; JV= -x 01 10, (4.4) 

Pio l) \ pot 1/ 
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The consistency condition (4.3) leads to the following discrete system of four equations: 

\i — v 



Pox - Pio 

qox - qxo 

Xoi — Xio 
^01 - ^10 



(1 +pnq) 2 



fj, — V 



i + pxxq + i>xxx)pxx 

1 +Pxxq + ipxu)q 



1 +Pxxq 
/j, — v 

l + Pxxq 



1> 



(4.5) 



If all odd variables vanish, this system of difference equations reduced to a familiar two- 
component system [1, 2, 3]: 



{fi - u)p n (fi - v)q 

Pox - Pio = , , — — ?oi - ?io " 



1 +Pxxq 



1 +Pxxq 



(4.6) 



(Pox, qox, 4>ox, *ox] 



rn 



-> k 



(PXO, qiO, IpXO, KXo) 



Figure 2: A staircase initial value problem for (4-5) 

One can easily solve (4.5) either with respect to (p w , q w , ip w , x 10 ) or with respect to (p i, 
qox, ''Pox, xqx)- 111 this case, one can pose an initial value problem for (4.5) with initial conditions 
on a staircase as it is shown on Figure 2. For a given set of initial data on the staircase, a 
solution of the difference system (4.5) can be found recursively. 

Similar to [33] one can define the Elimination map and express any variable on the lattice 
in terms of a finite subset of the initial set of variables on the staircase (Fig. 2). It is clear, that 
these expressions are rational functions of the even initial variables and multi-linear function 
of the odd ones. 

Case B: If we combine two Darboux transformations of types (3.14) and (3.8) from the 
previous section: 

((J> + Pioq ~ V^io + A ip q\ I F + A i/j q \ 

-xxo 1 , N = -*bi 1 ; (4.7) 

Pio l) V p m 0/ 



7 



then the compatibility (the Bianchi commutativity) condition will lead to the following quadri- 
lateral system: 



Poi = (A* - F + PioQ ~ *l>Kio)pu 
9io = (n~ F w + pugoi - ^01^11)9 
4>oi — V^lo = —(A* — ^10 +Pii?oi - ^Ol^ll)^ 
xbi - x 10 = (/i - F + pi g - ^xi )xn 



(4.8) 
(4.9) 
(4.10) 
(4.11) 



F{fi + Puqoi - i>aixu) = F 10 (fi + p 10 q - il>x 10 ) + pi q 10 -£>qi9oi - ^io><io + ■0oi><bi{4.12) 

and the condition: p i9 = 1 (3-9) which enables us to eliminate the variable p. One can solve 
(4.8) and (4.9) with respect to F and its shift F 10 : 

q 



T-! 9l0 . 

-r = A 4 1- 

9 



/ 77 9oi . 9oi , 

V^xio, F 10 = fi 1 Voi^u- 

9 9io 



(4.13) 



Then, using (4.13), one can eliminate F from (4.10) and (4.11). The compatibility condition 
S(F) = F 10 and (4.10)-(4.11)read: 



q q 9-1,1 9i,o , , c it n 

1 h V-i.iK ~ Vki-i + // - //i = 0, 

9-i,o 9i,-i 9 9 

^1,0 - V^oi = — 
9 



9io 



1,-1 



6o- 



(4.14) 
(4.15) 
(4.16) 



Here we have introduced a new variable £ = x 0; i. The equation (4.12) is satisfied due to (4.13) 
and (4.14)-(4.16). 

After setting q = e v , where v is an even variable, one can easily recognize (4.14) as a 
non-commutative extension of the fully discrete Toda chain: 



-.in,— i— " 



+ e 



v-i,o—v 



J-»1,0 



+ V>-i,i£ - V^i ,-i + A*i - A* = 0) 



(4.17) 
(4.18) 

^-6,-1 = 6^-^10. (4.19) 

In the special case when all anti-commuting variables vanish, it reduces to the discrete Toda 
lattice [7, 3]: 

(4.20) 



-."1.— l—v 



+ e 



u-i.o— v 



,1>-U1,0 



V—V—l 1 



+ /ii - /i = 0. 



On the lattice each equation (4.14) - (4.16) can be represented by a graph: 



•o 



Eq. (4.14): 



o o 



x x 
x< 

o 



(4.21) 



Eq. (4.15): 



Eq. (4.16): » 



x x 

X 



(4.22) 



(4.23) 
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where even variables are denoted by solid dots •, the variables ip - with empty circles o and 
the variables £ with crosses x. To each ratio in (4.14) there corresponds a solid line. 



(k, m) 




(k + 2p, m — 2p — 1) 



Figure 3: An initial value problem for (4-14)- The black staircase gives the initial values set 
for the even variables while the gray parallelogram gives the initial values for the odd variables 
(denoted by: o = i/j and x = £. The notation <S> means that we have both i/j and £ given) 

For the commutative Toda lattice (4.20) one can solve an initial value problem with initial 
data given on the staircase Wq: 

W = {(k + n,m-n),(k + n,m-n-l)\nE {0, ...,2p}}. (4.24) 

In the case of equations (4.14) one also needs to define boundary odd variables. Taking some 
p e N we define a parallelogram W with boundaries 

W\ = {(k + 2n, m - n), (k + 2n, m - n - p - 1) | n G {0, ...,p}} (4-25) 
W 2 = {{k,m-n)\ne {0, ...,p + 1}}, (4.26) 
W 3 = {{k + 2p,m - n - p) \ n e {0, ...,p + 1}}. (4.27) 

The set of boundary variables tp 1 ^ are defined on W\ U W 2 and the boundary variables £^ are 
defined on W± U W3. One can show, that the variables on the boundary W\ U W2 U W3 and 
inside W can be expressed as rational functions of the even variables given on the staircase Wq 
inscribed into the parallelogram W and multi-linear functions of the odd boundary variables. 

Indeed, the system (4.14) with such initial boundary conditions can be solved by a finite 
sequence of iterations. For the first iteration we set all odd variables inside the parallelogram 
to zero and find the first approximation of even variables for all points of W. Then, using the 
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boundary conditions for odd variables, one can solve equations (4.23) to update the values of 
odd variables inside W. Starting from these data, we repeat the sequence of iterations. This 
sequence will stabilise after a finite number of steps since the solution (<7& 1)mi , V^mu 6fei,mi)> 
(fci,mi) G W is a multi-linear function of the odd boundary data. 

5 Conclusions 

In the present paper we studied integrable difference equations associated with Grassmann 
extensions of the nonlinear Schrodinger equation. We constructed two elementary Darboux 
transformations. As a result, newGrassmann generalisations of the Toda lattice and the NLS 
dressing chain are obtained. We should mention here that the noncommutative NLS dressing 
chain we obtained is a nonlocal one - although the commutative reduction will lead to the 
standard NLS dressing chain [21]. The Toda lattice dressing chain obtained here differs from 
the known super-symmetric Toda chains [38, 13, 5, 18]. 

Furthermore, we obtained difference integrable systems as a compatibility (Bianchi com- 
mutativity) of these Darboux transformations. Such systems can be viewed as Grassmamm 
generalisations of the difference Toda and NLS equations. These lattice systems have Lax pairs 
provided by the set of two consistent Darboux transformations. The corresponding Backlund 
transformations represent symmetries of the discrete (difference) systems. 

For the two discrete systems obtained, we formulated the initial value problems: in the 
first case (when two Darboux transformations with NLS-type dressing chains are combined) 
the well-posed initial value problem is defined (for both even and odd variables) by giving an 
initial profile along an infinite staircase within the (k, m) -lattice (Fig. 2). In the second case, 
the well-posed initial boundary value problem is defined for a parallelogram (of any size). The 
initial even variables should be given on a staircase inscribed in the parallelogram while the 
odd variables should be given on the boundary of the parallelogram (Fig. 3). The Bosonic 
limit reconstructs the results of [21]. 

The osp(l|2)-invariant supersymmetric NLS model of Kulish [ ] can be obtained from 
(2.3) by imposing the reduction p = q' and %fj = ( = <ft = xL Our Darboux transformations 
is not applicable to the Kulish system, because they do not respect the reduction to osp(l\2) 
super algebra. 

The results obtained here can be developed in several directions: 1) to study the correspond- 
ing Yang-Baxter maps [20, 39]; 2) To derive the recursion operators and study the associates 
multi-Hamiltonian structures; 3) This can be generalised to other integrable hierarchies. 
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